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Using a sample of 14 million ip{2S) events accumulated with the BES II detector, evidence for 
/o(980)/o(980) production in XcO decays is obtained for the first time; the branching ratio is deter- 
mined to be Bixco -» /o(980)/o(980) ^ Tr+Tr-Tr+Tr") = (7.6 ± 1.9 (stat) ± 1.6 (syst)) x 10"''. The 
significance of the /o(980) signal is about 4.6a. 

PACS numbers: 13.25.Gv, 14.40.Cs 

I. INTRODUCTION 

After thirty years of controversy, the nature of the /o(980) is still not settled [1]. It has been described as a 
conventional qq meson [2] , a "unitarized remnant" of a qq state [3] , a KK molecule [4] , a multiquark state [5] , or a 



glucball [6]. Because of its close proximity to the KK threshold and its propensity to decay to KK^ it is difficult 
to quantify even the mass and width of the /o(980). To be explicit, the state with a mass of 980 ± 10 MeV and a 
width somewhere between 40 and 100 MeV [7] straddles the KK threshold at 990 MeV. Many arguments favoring or 
disfavoring the above assignments depend on the width or pole position of the /o(980). 

A novel measurement to elucidate the nature of the /o(980) was suggested by Refs. [8,9]. By determining the 
radiative decay rate for cf) -^ /o(980)7, one can infer the ss content of the /o wave function since the rate is proportional 
to the overlap with the </), a well known ss state. The results from CMD2, SND, and KLOE [10] reveal a much higher 
branching ratio for radiative (j) —^ 7/0 decay than that expected for the qq meson or KK molecule interpretations. It 
seems that these data add weight to the idea that the /o(980) is a compact qqqq state with an extended meson- meson 
cloud 'molecular' tail [11]. However, at present the interpretation about the nature of the /o(980) is still open [7], 
and more experimental results are needed to clarify it. 

In this paper, we report on the analysis of 7r+7r^7r+7r^ final states from Xco decays using a sample of 14 million 
tp{2S) events accumulated with the BESII detector. Evidence for /o(980)/o(980) production from XcO decays is 
obtained for the first time. 

II. BES DETECTOR 

BES II is a large solid-angle magnetic spectrometer that is described in detail in Ref. [12] . Charged particle momenta 
are determined with a resolution oiap/p = 1.78%-\/l + p^ {p in GeV/c) in a 40-layer cylindrical drift chamber (MDC). 
Particle identification is accomplished by specific ionization (dE/dx) measurements in the drift chamber and time-of- 
flight (TOF) measurements in a barrel-like array of 48 scintillation counters. The dE/dx resolution is (idE/dx = 8.0%; 
the TOF resolution is (Ttof = 180 ps for Bhabha events. Outside of the timc-of-flight counters is a 12-radiation- length 
barrel shower counter (BSC) comprised of gas proportional tubes interleaved with lead sheets. The BSC measures the 
energies of photons with a resolution of oe/E ~ 2\%/^/E {E in GcV). Outside the solenoidal coil, which provides a 
0.4 Tesla magnetic field over the tracking volume, is an iron flux return that is instrumented with three double layers 
of counters that are used to identify muons. 

In this analysis, a GEANT3 based Monte Carlo simulation package (SIMBES) with detailed consideration of detector 
performance (such as dead electronic channels) is used. The consistency between data and Monte Carlo has been 
checked in many high purity physics channels, and the agreement is quite reasonable. 

III. EVENT SELECTION 

The selection criteria used in this analysis are similar to those of Ref. [13]; the main difference between them is 
that no particle identification is imposed here in order to increase the selection efficiency. 

A. Photon Identification 

A neutral cluster is considered to be a photon candidate when the angle between the nearest charged track and 
the cluster is greater than 15°, the first hit is in the beginning 6 radiation lengths, and the difference between the 
angle of the cluster development direction in the BSC and the photon emission direction is less than 30°. The photon 
candidate with the largest energy deposit in the BSC is treated as the photon radiated from the ip{2S) and used in a 
four-constraint kinematic fit to the hypothesis ip{2S) -^ 77r'^7r~7r+7r~. 

B. Charged Particle Identification 

Each charged track is required to be well fit to a three-dimensional helix using the MDC information, be in the 
polar angle region | cos9mdc\ < 0.80, and have the point of closest approach of the track to the beam axis be within 
2 cm of the beam and within 20 cm of the center of the interaction region along the beam line. 

C. Event Selection Criteria 

The candidate events are required to satisfy the following selection criteria: 



1. The number of charged tracks is required to be four with net charge zero. 

2. The sum of the momenta of the lowest momentum 7r+ and n^ tracks is required to be greater than 650 MeV; 
this removes contamination from ip{2S) —> tt+tt^ J/i/; and some p'^mr events. 

3. The x^ probabihty for the four-constraint kinematic fit to the decay hypothesis ip{2S) -^ 77r+7r^7r+7r^ is greater 
than 0.01. 

The invariant mass distribution for the n'^n^Tr^T:^ events that survive all the selection requirements is shown in 
Fig. 1. There are clear peaks corresponding to the Xcj states. The highest mass peak corresponds to charged track 
final states that are kinematically fit with an unassociated photon. 
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FIG. 1: The 7r"'"7r~7r"'"7r~ invariant mass spectrum. The shadow histogram shows the spectrum for Monte Carlo simulated 
background events. 

The distribution of background events in the Att mass spectrum, determined by Monte Carlo simulation normalized 
using PDG2004 branching ratios [7], is also shown in Fig. 1. The distribution is flat in the mass range of the 
Xcj states, and the background events come mainly from il^{2S) -^ 7r'^7r+7r^7r+7r^. The highest mass peak is from 
'0(25') — > TT+TT^TT+Tr" and '0(25') -^ TT^Tr~ K'^ K~ events combined with an unassociated photon. The background is 
very low compared with the strong Xcj peaks, and its effect will not be considered in the following analysis. 

In this analysis, no particle identification is imposed on charged tracks in order to increase the detection efficiency. 
It is not necessary to distinguish pious from kaons or protons in this channel because the background is not serious, 
as shown in Fig. 1, and the contamination from events with kaons or protons is rejected effectively by the kinematic 
fit. The statistics of present data sample is about 20% higher than one using particle identification. 



IV. ANALYSIS RESULTS 



A. /o (980) /o (980) signal 



Figure 2 shows scatter plots of tt+tt^ versus tt+tt" invariant mass [14] for events with 7r+7r^7r+7r^ mass between 
3.30 and 3.48 GeV and between 3.53 and 3.60 GeV, and the corresponding projections are shown in Fig. 3 (two 
entries per event). The clusters of events in the lower left-hand corners of Figs. 2(a) and 2(b) indicate the presence of 
a KgKg signal under both the Xco and Xc2 peaks. A clear /o(980)/o(980) signal can be seen in Fig. 2(a). There are 
some hints of pV° and /o(1370)/o(1370) (or /2(1270)/2(1270)) signals in Fig. 2(a) and p°p° but no /o(980)/o(980) 
events in Fig. 2(b). In this paper, we study the /o(980)/o(980) in XcO decays. 

For the events in XcO mass region (from 3.30 to 3.48 GeV) and after requiring that the mass of one of the tt+tt^ 
pairs lies between 0.88 and 1.04 GeV, the mass distribution of the other tt+tt" pair is shown in Fig. 4 (two entries 
per event); there is a strong /o(980) signal, and its line shape is similar to other experiments [7]. From a Monte 
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FIG. 2: Scatter plots of tt'^tt versus tt^tt invariant mass for selected j-k'^tt ■k'^-k events with tt'^tt ■k'^-k mass in (a) the 
Xco and (b) the Xc2 mass regions. 
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FIG. 3; Projections of tt+tt^ invariant mass under the (a) XcO and (b) Xc2 peaks (two entries per event). 
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FIG. 4; Plot of -K^-K mass recoiling against the /o(980) (0.88 GoV < ra^^^- < 1.04 GcV) for events in the XcO mass region 
(two entries per event), where the dashed line histogram indicates a rough estimation of background determined from sidebands. 
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FIG. 5: Scatter plot of vr tt versus n 
of signal and background regions. 



TT invariant mass in the /o(980) region for XcO candidate events, showing the definition 



Carlo simulation, the background in the /o(980) region is mainly from processes such as ip{2S) —^ "fXcO, XcO ~* 
ai(1260)7r, ai(1260) ^ pvr. 

Note that the background estimation in Fig. 4 using sidebands (0.76 - 0.84 GeV and 1.08 - 1.16 GeV) is rough. 
In this paper, the number of /o (980) /o (980) events and the corresponding background are estimated from the scatter 
plot of TT+TT" versus tt+tt" invariant masses, as shown in Fig. 5. This method gives more accurate determinations of 
the /o(980) signal and background. The signal region is shown in Fig. 5 as a circle centered at (0.960,0.960) GeV 
and with a radius of 80 MeV, and the background is estimated from the events between two circles with radii of 120 
MeV and 160 MeV. There are 65 and 51 events in the signal and background regions, respectively. So the number 
of /o(980)/o(980) events is estimated to be 65 - 51/1.75 = 35.9 ± 9.0, where 1.75 is the normalization factor - the 
ratio of the area of background region to that of the signal region. The tt+tt" mass range we adopt is shifted from 
the /o(980) central mass value of 980 MeV because of the asymmetric character of its mass spectrum. We obtain the 
signal significance of the /o(980)/o(980) of 4.6cr using the method described in Ref. [15]. 



B. Monte Carlo simulation 



A Monte Carlo simulation is used to determine the detection efficiency. The angular distribution of the emitted 
photon in the process ip{2S) -^ 7XcO is taken into account [16]. The /o(980) is generated with the usual Flatte formula 
[17, 18]: 



./ = 



1 



A/2 



iM{giPi + g2P2) ' 



where pi2 are the phase space factors for the tttt and KK channels, Pi{s) = -y/l — mf/As, and g\_2 are squares of 
coupling constants to the two channels. For the KK channel, m? is taken as the average of the K'^ and K^ masses, 
and the algebraic expression for p2 is extended analytically below the KK threshold. In the simulation, the parameters 
used are those of Ref. [18]: M = 0.9535, .gi = 0.1108, 32 = 0.4229 GeV. 



C. Branching fraction results 



The efficiency is determined using 100,000 Monte Carlo simulated events that are passed through the same selection 
as the data events; the efficiency is estimated to be e = (3.92 ± 0.07)%, where the error is the statistical error of 
the Monte Carlo sample. Note that for this estimation, the events in the background region are subtracted from the 
events in the signal region, similar to the treatment of data. 



Using numbers from above, the branching ratio of ip{2S) -^ 7XcO: Xco -^ /o(980)/o(980) -^ tt+tt vr+vr is 
B{i^{2S) -> 7XcO -> 7/o(980)/o(980) ^ 77r+7r-^+^-) = (6.54 ± 1.64) x l{r\ 
where the error is statistical. 
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FIG. 6: A fit to the -k^t^ mass recoiling against the /o(980) (0.88 GcV < m^+^- < 1.04 GeV) for events in the XcO mass 
region (two entries per event). 

As a cross-check, a fit to the projected tt+tt^ mass distribution between 0.49 and 1.81 GeV with the Flatte 
parameters of the /o(980) fixed to the solution of Ref. [18] plus a polynomial background yields 84.8 ± 18.4 events for 
the /o(980) signal, as shown in Fig. 6 (two entries per event). From the results of this fit, we determine the branching 
ratio 



3^23) ^ 7Xco ^ 7/o(980)/o(980) 



77r^7r 



TT+TT-) = (4.54 ±0.98) X 10" 



where the detection efficiency is 13.34% (two entries per event). Because of the low statistics and the relatively high 
backgrounds, as well as the lack of information from the coupled KK channel, a free fit to the parameters M, gi, and 
g2 is difficult. The branching ratio values determined by these two methods agree to about one sigma. 

D. Systematic errors 

The systematic errors in the branching ratio measurement associated with the efficiency are determined by com- 
paring 4'{2S) data and Monte Carlo simulation for very clean decay channels, such as ip{2S) -^ tt+tt^ J/i/i, which 
allows the determination of systematic errors associated with the MDC tracking, kinematic fitting, and the photon 
identification [19]. Other sources of systematic error come from the uncertainties of the number of tpi'^S) events [20], 
the parameters of the /o(980), the definition of background region, the XcO and /o(98G) mass resolutions, etc. 



Parameters of the /o (980) 



The parameters of the /o(980) are still uncertain, and different descriptions of the /o(980) in the simulation result in 
different efficiencies. Besides the solution of Ref. [18], we also consider the measurements of some recent experiments 
such as E791, GAMS and WA102 [21-23], where a Breit-Wigner description with the width varying from 44 to 80 
MeV was used for the /o(980). We determine the change both by using the solutions of Refs. [21-23] and by varying 
gi in Ref. [18] from 0.1108 GeV to 0.090 GeV and 0.130 GeV while keeping the ratio 52/31 fixed. The largest change 
is about 16%, which is used for the systematic error due to this uncertainty. 



2. Different background regions 

In this paper, we estimate the background using the region between circles with radii 120 MeV and 160 MeV about 
(0.960,0.960) GeV, as shown in Fig. 5. We test two other different background definitions by changing the radii of 
the two circles to 100 and 150 MeV and 120 and 180 MeV. The biggest change is about 5%, which is taken as the 
systematic error. 

3. Xco <J-nd /o(980) mass resolutions 

Differences between data and Monte Carlo mass resolutions for the XcO and /o(980) also cause systematic uncertain- 
tics in the determination of the branching ratio of XcO ~^ /o(980)/o(980). From a study, we find that the difference for 
the XcO is about 1 MeV, so we change the window of XcO to [3.300 + 0.005, 3.480 - 0.005] GeV and [3.300 - 0.005, 3.480 
+ 0.005] GeV, and estimate the effect on the branching ratio. Such changes result in less than a 1% variation in the 
efficiency, and the effect of the difference in the mass resolutions of the /o(980) is even smaller. By varying the width 
of XcO by Icr of its error, 0.8 MeV, there is almost no change on the detection efhciency. We include a 1% systematic 
error for the sum of these uncertainties. 

Table I lists the systematic errors from all sources, and adding them in quadrature, the total systematic error, 20%, 
is obtained. The resulting branching ratio is 

B{i;{2S) -^ 7XcO ^ 7./o(980)/o(980) ^ jtt+tt-tt+tt-) ^ (6.5 ± 1.6 ± 1.3) x 10-^ 

and finally using the PDG2004 average value and error for B{ip{2S) -^ 7Xco) [7], we obtain 

^(Xco -^ /o(980)/o(980) -^ 7r+7r"7r+7r") = (7.6 ± 1.9 (stat) ± 1.6 (syst)) x lO"''. 



TABLE I; Summary of systematic errors in the branching ratio calculation of B{tp{2S) — > 7XcO -^ 7/o(980)/o(980) 

77r TT^TT TT^). 



Source 


Relative systematic error 


MDC tracking 


8% 


Kinematic fit 


6% 


Photon ID efficiency 


2% 


'ip{2S) number 


4% 


Efficiency estimation 


16% 


Definition of background 


5% 


Mass resolutions 


1% 


Total 


20% 



V. SUMMARY 

Evidence for /o(980)/o(980) production from XcO decays is obtained for the first time with a significance of about 
4.6(7 , and the branching ratio is determined to be B{xcO -^ /o(980)/o(980) -^ 7r+7r^7r+7r^) = (7.6 ± 1.9 (stat) ± 
1.6 (syst)) X 10"**. This may help in understanding the nature of /o(980). 
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